Non-equilibrium Luttinger liquid: Zero-bias anomaly and dephasing 
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A one- dimensional system of interacting electrons out of equilibrium is studied in the framework 
of the Luttinger liquid model. We analyze several setups and develop a theory of tunneling into such 
systems. A remarkable property of the problem is the absence of relaxation in energy distribution 
functions of left- and right-movers, yet the presence of the finite dephasing rate due to electron- 
electron scattering, which smears zero-bias-anomaly singularities in the tunneling density of states. 



PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 

The interest in one-dimensional (ID) interacting elec- 
tron systems is due to their fascinating physical prop- 
erties and potential applications in nanoelectronics. A 
variety of experimental realizations of quantum wires 
includes carbon nanotubes, quantum Hall edges, semi- 
conductor structures, polymer fibers, and metallic 
nanowires. At equilibrium, the physics of ID electrons 
has been thoroughly explored. The main peculiarity of 
the system is the formation of a strongly correlated state, 
Luttinger liquid (LL), commonly described in terms of 
collective bosonic excitations [l[ • A hallmark of LL cor- 
relations is a strong, power-law suppression of tunneling 
current at low bias — zero-bias anomaly (ZBA) [2|]. 

On the other hand, little is known about LL away from 
equilibrium. Theoretical efforts so far focused on non- 
linear transport through a single impurity [3, 0] . In this 
work we consider the problem of tunneling into a LL out 
of equilibrium. It is important to emphasize a remark- 
able property of the LL: the absence of relaxation towards 
equilibrium. Once non-equilibrium energy distributions 
of left- and right-movers are created, they propagate 
along the wire arbitrarily long without thermalization 
(within the clean LL model, i.e. when back-scattering 
by interaction or impurities, spectral nonlinearity, and 
momentum dependence of interaction are neglected) [B[ . 
Below we address such a situation and analyze the tun- 
neling density of states (TDOS) of a non-equilibrium LL 
wire. We show that while the energy relaxation rate is 
zero, the ZBA in TDOS is smeared out by dephasing 
processes yielding a finite quasiparticle life time. 

We begin by discussing possible experimental realiza- 
tions of a non-equilibrium LL. Throughout most of the 
paper we assume that the interaction strength interpo- 
lates adiabatically between its value in the LL (central 
part of the wire where the measurements are performed) 
and zero near the electrodes. The assumption of adia- 
baticity is not entirely innocent; at the end, we briefly 
discuss the case of sharp switching of interaction and en- 
suing modifications. 

The simplest setup is shown in Fig. QJi. A long clean 



LL is adiabatically coupled to two electrodes with differ- 
ent potentials, fir, — (J-r = eV and different temperatures 
T v (where rj = L,R stands for left- and right-movers) 0]. 
A particularly interesting situation arises when one of 
temperatures is much larger than the other, e.g., T L = 
and Tr finite. Then the ZBA at fir,, where the distri- 
bution function has a sharp step, is broadened solely by 
dephasing originating from electron-electron scattering. 




FIG. 1: Schematic view of setups for measurement of the 
TDOS of a LL out of equilibrium (see text for details). 

A more complex situation arises if electrons coming 
from reservoirs with different potentials mix by impurity 
scattering. Two different realizations of such devices are 
shown in Figs. [T]3,c. In the first case, Fig. QJ>, the mix- 
ture of left and right movers coming from reservoirs with 
Mi is caused by impurities which are located in the 

non-interacting part of the wires. In the second setup, 
Fig- Ht) the LL wire is attached to two metallic wires 
which arc themselves biased. We assume that these elec- 
trodes are diffusive but sufficiently short, so that inter- 
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electrode energy equilibration can be neglected. As a re- 
sult, a double-step energy distribution is formed in each 
electrode @. The left- and right- movers in the LL wire 
"inherit" these non-equilibrium distributions emanating 
from the respective electrodes. The existence of multiple 
Fermi edges in the distribution functions "injected" from 
the electrodes renders the behavior of the TDOS highly 
non-trivial, since the ZBA is expected to be broadened 
by electron-electron scattering processes governing the 
dephasing rate r^. 

The question of non-equilibrium ZBA broadening in- 
duced by electron-electron scattering is particularly in- 
triguing in the case of a ID system. First, energy relax- 
ation is absent in that case. Second, there are two quali- 
tatively different predictions concerning dephasing in the 
context of weak localization and Aharonov-Bohm oscil- 
lations: while the weak-localization dephasing rate van- 
ishes in the limit of vanishing disorder [9( , the Aharonov- 
Bohm dephasing rate is finite in a clean LL 0, [l(| • It 
is thus very interesting to see how dephasing processes 
manifest themselves in the broadening of ZBA. From the 
technical point of view, the challenge is to develop meth- 
ods to treat LL away from equilibrium. 

Within the LL model, the electron field is decoupled 
into the sum of left- and right-moving terms, ip{x,t) = 
ip R (x, t)e tpFX + ip L (x, t)e~ ipFX . The Hamiltonian reads 



H 



where Vq is the bare electron-electron interaction. We 
will neglect the interaction between the tip and the wire, 
assume that the tunneling current is weak and that elec- 
trons in the tip are at equilibrium at a negligibly low 
temperature. We further assume that the TDOS of the 
tip, i>ti can be considered as constant. Then, dependence 
of the differential tunneling conductance on the voltage 
Vt at the tip is controlled by the TDOS of the wire 



dl/dVt oc \T\ 2 v t v{eV t ) 



(1) 



where T is a tunneling matrix element. We thus focus 
on the TDOS v(e) of a non-equilibrium LL. 

In order to find the TDOS v{e) = v R (e) + u L (e) 
at a point x one needs to calculate the retarded 
("r") and advanced ("a") single particle Green func- 
tion, v n (e) = ^[G r v (e,x,x) - G«(e, x, x)]. To do it 
away from equilibrium we employ the Keldysh formal- 
ism. The Keldysh Green function reads G v (x, t; x' , t') = 
— i (Tk^ v (x, t)ij)Ux', t')), where the time ordering is 
along the Keldysh contour. We proceed by following 
the lines of functional bosonization approach [11 , 12, 13] • 
While fully equivalent to the conventional bosonization 
technique for the case of clean equilibrium LL, this 
method is advantageous in more complicated situations. 
In particular, the efficiency of the functional bosonization 
for the analysis of transport and interference phenom- 
ena in a disordered LL was recently demonstrated Q. 



The key features of the functional bosonization are: (i) 
it retains explicitly both fermionic and bosonic degrees 
of freedom; (ii) for Vq = it straightforwardly reduces 
to the model of non-interacting fermions. This makes 
the functional bosonization an appropriate tool for the 
development of the theory of a non-equilibrium LL. 

Decoupling the interaction by Hubbard-Stratonovich 
transformation via a bosonic field 6, we obtain the action 



(2) 



where d R L = 9* ± vd x . It is convenient to perform a 
rotation in Keldysh space, thus decomposing fields into 
classical and quantum components, ip = (i/'+±V'-)/v / 2 
(where + and — label the fields on two branches of the 
contour) and analogously for </>. We further introduce 
vector notations by combining cj> and <f> in a 2-vector (j). 

The Green function of interacting electrons can be pre- 
sented in the form G v (t — t',x — x') = (G v< j,(x, t; x' , t')), 
where G n< j> is the Green function of non-interacting 
fermions in an external field (f>, the averaging goes with 
the weight Z^e - ^ v ° and is a sum of vacuum 

diagrams (fermionic loops) in the field 4>. The special 
feature of ID geometry is that the coupling between 
the fermionic and bosonic fields can be eliminated by a 
gauge transformation, ip v (x, t) — > ip v (x,t)e* e i( x,t ', with 
@« = + o'iO, if we require 



id n 9. n 



(3) 



As a result, G^ can be cast in the form 

G^ix, t; x',t') = e^^Grjoix, x';t- t'y-' 1 ®^ , 



""' 1 Gf l0 



(4) 



Here G^o is a Green function of free fermions, with 
the Keldysh component G*f (e) = [1 — 2n r ,(c)][GZ Q (e) — 
G^o(e)], and n.^(e) is fermionic distribution function. 

To proceed further, we use the random-phase ap- 
proximation (RPA), within which Z^ is Gaussian, 
logZfi = — g <p II<p. This is an exact relation at equi- 
librium 1J], which is crucial for the LL being an ex- 
actly solvable problem. It remains exact for the non- 
equilibrium setup of Fig. [T^,, where both distributions 
n v are of Fermi-Dirac form. On the other hand, RPA 
becomes an approximation for more general non-equilib- 
rium situations (cf. Figs. []J>,c); we discuss its status and 
possibility of an exact solution at the end. 

The polarization operator of free fermions is given by 
LT = U R + U L , with 
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2tt uj + =p vpq 

(n;-ir^>). 
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where lj± — u> ± iS. The function 



B v Jlj) = - / den v (e)[2-n v (e-cj)-n v (e + uj)}, (6) 

w J-oo 

is related to the distribution function N^(u>) of electron- 
hole excitations moving with velocity v in direction r], 
B"{lu) = l + 2AP(w). At equilibrium, B*(w) = B cq {uj) = 
l + 2N eq (ui), where N cq (ui) is the Bose distribution. Since 
fcrmions are free after the gauge transformation, no re- 
laxation of the distribution functions n v (e) takes place. 

Performing the averaging over 0, we express the TDOS 
in terms of the correlation function of the gauge fields 



^ = 1 + 2, 
vq 



dtn^e^p (-1^8^(1^). (7) 



Here vq is the bare (non-interacting) density of states, 
n^it) is the Fourier transform of n^(e), 



C\t) = / (du,)(dq)(l - cosM) (69)W Q e-M/ A 



1%' =2 (dw)(dq) S m(u;t)(66) ™e 



(8) 



and A is an ultraviolet cutoff. To calculate Igg and I e g 
one needs to resolve Eq. ([3]) and to express the gauge 
field 6 in terms of the Hubbard-Stratonovich filed <b, 



(9) 



where GrjO is the Green function of free bosons with the 
Keldysh component Q^ = (£^ — G^ a )B" r The latter 
serves to reproduce correctly the distribution function 
in the (66) correlation function, ensuring, in partic- 
ular, the fluctuation-dissipation theorem at equilibrium. 

The correlation functions of 6 fields can be now readily 
found. For the (60) component (which is independent of 
the distribution functions) we get 



o)i R q L) = = 



iV 



uj ± vq 



2(tu + ^fvq) 2 2 (u> + =p vq)(uj\ — u 2 q 2 ^ 

where u = v(l + Vq/ttv) 1 / 2 is the sound velocity. 
This yields for the g-integrated propagator [which enters 
Eq. ©] J(dq)(00) { ^ = ^ 7 /2w, where 7 = (K - 1) 2 /2K 
and K — v/u = (1 + Vq/ttv)^ 1 / 2 is the conventional di- 
mensionless parameter characterizing the LL interaction 
strength. Evaluation of the u> integral leads to 



Af(t) = 7arctan(l + t 2 A 2 ) 



(10) 



At equilibrium / (dq)(00) { $ = (^7 /uj)B cq (uj). For T = 
this yields 



4^) = |log(l + * 2 A 2 ) 



(11) 



Substituting Eqs. (jTTJ) , (fTT)|) into Eq.((7]), we reproduce 
the famous power-law behavior of TDOS, 



^cq(e) ~ vq (e/A)' 



(12) 



At finite T the long-time behavior of iil) is modified, 



I$(t)~^Tt = t/2 T<l> . 



(13) 



The ZB A dephasing rate is thus 1 /t$ = 2irjT and agrees 
with the one relevant to Aharonov-Bohm effect 0, [lj| • 
While I/t^ contributes to smearing of ^ eq (e), for 7 < 1 
this is not particularly important, as the singularity is 
anyway smeared on the scale of T due to the distribution 
function n^t) in Eq. l[7]). 

We turn now to the non-equilibrium situation. It is 
important to emphasize that the functional bosonization 
procedure can be also applied to a (clean) LL with spa- 
tially dependent interaction constant K(x) (such a model 
was considered in [HI]). Since the interaction can still be 
gauged out, the conclusion about independence of the 
fermionic distribution n v (e) on x retains its validity. On 
the other hand, the boson sector will be, most generally, 
characterized by different distribution functions B™(u>), 
B^(u>) corresponding to the modes propagating with ve- 
locities u and v [coexistence of such modes is clear from 
the correlation function (|10p]. While the u-excitations 
are conventional plasmons, the w-mode describes bare 
electron-hole pairs (thus moving with velocity v of non- 
interacting fermions). The corresponding formalism for 
higher-dimensional diffusive system was developed in (l6j 
(where the analog of the w-mode was termed "ghosts" ) ; 
for an analysis of energy relaxation in disordered LL in 
this framework, see Ref. \vA. We find 



f (dq)(60) i 



uj\q 



■(7 + 2)^0 



and analogously for 

While the distribution B~l of bare electron-hole pairs 
is determined by the fermion distribution function n, q , 
see Eq. ([6]), the plasmon distribution B™ may differ, de- 
pending on the setup. We consider first the "adiabatic" 
situation when K (x) changes slowly on the scale 1$ deter- 
mined below. Then plasmons with relevant wave vectors 
are not backscattered by modulation of the interaction 
K(x), so that m = B v v = B n is given by Eq. ©. In this 
situation Eq. (fQ|) reduces to 



(dq)(00)* = ih[Br{u) + B L (ou)]/2u . (14) 



For the setup of Fig. [Qa, Eq. dHJ) yields I { e f = nj(T L + 
T R )t/2. Therefore, the ZBA dips in TDOS of both chiral 
sectors (separated by fi^ — Mi? : = eV) get broadened (in 
addition to the thermal smearing) by the dephasing rate 



1/% = K7( T L + T R ) . 



(15) 
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FIG. 2: ZBA in LL, setups b,c with a = §, A = 1, T,, = 0, 
eV = 0.1, for 7 =0.1, 0.2, 0.25 (from top to bottom). 



If one of Tft t L is much lower than the other, so that e.g. 
Tl <C the broadening of the corresponding ZBA 

minimum (near fir,) is determined by the dephasing. 

For the "fully non-equilibrium" setups of Fig. [Tj3,c the 
fermionic distributions n v (e) are of a double-step form, 



n v (e) = a n n (e-) + (1 - a 7) )n (e + ) , 



(16) 



where no(e) is the zero-T Fermi distribution (we assume 
T n < eV), 0<a v <l, and e± = e ± U/2 . For the distri- 
bution (|16j) . one finds the bosonic distribution function 

ujB^(uj) = [a* + (1 - a v ) 2 ]ujB cq (uj) + a v (l - a v ) 
x [(ui + eV)B cq (uj + eV) + (u - eV)B cq {u - eV)] . 

This yields the TDOS (see Fig. |2J) 

v v (e) ~ a^ eq (max{e_, 1/2-r^}) 



with the non-equilibrium ZBA dephasing rate 
l l T l = c v irjeV 



(18) 



and the numerical prefactor c v — or(1 — or)+<Z£(1 — a^). 

If iiC(^) varies fast, the plasmon distribution be- 
comes spatially dependent, while remains unchanged. 
In the limit when the interaction is turned on as a sharp 
(on the scale 1$) step, the plasmons are scattered accord- 
ingly to the Fresnel law [15j], with a reflection coefficient 
R=(l — K) 2 /(l + K) 2 . This yields a boundary condition 
for the distributions B™; e.g. for a sharp change on the 
r.h.s. of the wire B£(u;)=(l - R)B v l {lo) + RB^uS), and 
similarly for the left boundary. In this way it is easy to 
treat the situation with sharp switching of the interaction 
on one or both sides [2(J. The results (fl7|) . (fl"8|) retain 
their validity but with modified numerical prefactors c v . 

To summarize, we have analyzed several setups in 
which non-equilibrium LL can be observed. Using the 
functional bosonization formalism, we have developed a 



theory of tunneling into such systems. While energy re- 
laxation is absent, the dips of split ZBA are broadened 
by dephasing, Eqs. (fl~5|) . (fl~8|) . yielding a finite quasipar- 
ticle life time. We reiterate that RPA is exact for the 
setup of Fig. [Iji, but not for the setups of Figs. [Tj3,c. In 
particular, the RPA value of the prefactor c v in Eq. (fT8"|) 
for I/t^ should be considered as an approximation. One 
can expect that RPA becomes controllable for weak in- 
teraction, 7 -C 1. Also, it remains to be seen whether ex- 
act results can be obtained for a generic non-equilibrium 
setup. The key observation is that the sum of vacuum 
diagrams (log Z$) can be cast in the form analogous to 
the gen erating function in the counting statistics problem 
Til . l^ | . Work in these directions is underway. 
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